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Abstract 

We study the effect of a close encounter of a passing star on the shape of the inner 
Oort Cloud, using the impulse approximation. The deviation of the perturbed Oort 
Cloud from sphericity adds angular fluctuations to the brightness of the Cosmic Mi- 
crowave Background (CMB) due to thermal emission by the comets. An encounter 
with a solar an impact parameter of 2100 AU, as expected based on 

the abundance and velocity dispersion of stars in the solar neighborhood, leads to 
a quadrupole moment C2 = 4.5(3.5) x 10-l^ 6.7(1.1) x lO^^^ 1.1(0.11) x lO"*^ at 
u = 30, 353, 545 GHz, respectively in intensity and (temperature) fluctuations. We 
also quantify the quadrupole spectral distortions produced by the Scattered Disc, 
which will exist irregardless of any perturbation and the subsequent shape of the 
Oort Cloud. For comparison, the temperature fluctuation quadrupole moment pre- 
dicted by the current cosmological model is C2 = 1-76 x lO"^'^, which corresponds 
to fluctuation in the CMB intensity of C2 = 2.9 x 10-i°,6.8 x 10-^,1.6 x 10"* 
at v = 30, 353, 545 GHz. Finally, we discuss how a measurement of the anisotropic 
spectral distortions could be used to constrain the trajectory of the closest stellar 
fly-by. 
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1 Introduction 



At the outer edge o f the Solar System exists a swar m of small icy rocks known 
as the Oort Cloud (Oort, 19501 : Pones et al. . 2004 ). Comets with long orbital 



Email address: babich@tapir.caltech.edu (Daniel Babich). 



Preprint submitted to Elsevier 



23 July 2008 



periods, P > 200 yr, are ejected from this cloud into the inner Solar System by 
the perturbations from passing stars, giant molecular clouds or tides from the 
Galactic disc. The perturbed comets are moved onto nearly-parabolic orbits 
which approach the Sun, develop comae and may be detected. Gravitational 
perturbations by the planets can strongly alter the comets' orbits as they enter 
the inner Solar System, and so typically they either get ejected from the Solar 
System or left strongly bound to the Sun. It is very unlikely that a comet will 
return again on a nearly parabolic orbit. The observe d con t inual flux of long 
period comets on nearly parabolic orbits originally led lOort fll950f ) to propose 



that all long period comets were entering the inner Solar System for the first 
time and therefore a resevoir of comets must exists in the outermost region of 
the Solar System. 

There is considerable interest in studying the Oort Cloud since it is a rem- 
nant of the formation epoch of the Solar System. This is true for both the 
individual objects, whose chemical composition may reveal information about 
the compositi on and thermal state of the outer regions of the proto-planetary 



disc (but see iMumma et al.l (Il993l ) for a comprehensive review of processes 
that may cause significant evolution in these properties over the past 4 Gyr) 
as well as the dynamical structure of the Oort Cloud as a whole. The cloud 
structure, which is largely determined by its formation scenario, may retain 
information about the masses of the giant planets, the rate at which passing 
stars perturbed the comets' orbits and the surface density of any gas that 
was potentially still within the Solar System. This information could provide 
important clues about the formation processes of the outer planets. 

The Oort Cloud is cruedly divided into an inner and an outer region. This 
division stems from the fact that the probability for ejection into the inner 
Solar System by externa l perturbations is a stro ng function of semi-major 



axis (IHeisler &: Tremaind . Il986l : iDones et al.l . 120041 ) . All of the observed long- 



period comets are believed to have originated from the outer Oort Cloud. The 
perturbations that eject comets into the inner Solar System also isotropize 
the comets' orbits and eliminates any information about the formation epoch 
in the structure of the outer Oort Cloud. Therefore the reduced sensitivity 
of the comets in the inner Oort Cloud to these external perturbers is both 
good and bad. On the one hand it leads us to suspect that the structure of 
the inner Oort Cloud should retain an imprint of its formation as well as 
any rare, strong gravitational scattering event that might have occured over 
the past 4.5 billion years. On the other hand, it makes studying the inner 
Oort Cloud difficult. It is impossible to directly observe the comets in the 
inner Oort Cloud through their refiected sunlight because they are too far. 
Extremely large ob jects like Sedna can be optically detected when they are 



close to perihelion (IBrown et al.l . 120041 ). however it will be difficult to detect 



the potentially large number of small objects {R < 1 km) in this region. One 
novel technique involves using the microwave emission from these bodies to 
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constrain, or possib ly detect, the bulk properties of objects in this region 
flBabich et al1 . l2007h . 



At a distance D from the sun an Oort Cloud object will be heated by the Sun 
to an equilibrium temperature, 



Oort 
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1000 AU 



D 



(2) 



here we adopted an albedo of A = 4% (jjewitt fc Lud . Il998l ) and an emissiv- 
ity of e ~ 1.0. We adopt a model such that the emissivity of an Oort Cloud 
object is approximately unity when its size is larger than the wavelength of 
the relevant radiation and zero when the wavelength is larger. Also, a is the 
Stefan-Boltzmann constant and Lq = 3.83 x 10^^ ergs is the Solar luminos- 
ity. The equilibrium temperature is higher than t he Cosmic Mi c rowav e Back- 
ground (CMB) temperature of Tcmb = 2.725 K dMather et all . Il999h . When 
the comets are farther away from the Sun a variety of physical pro cesses help 
maintain temperatures greater than 5 — 6 K (iMumma et al.l . Il993l ). The Oort 
Cloud objects extinguish some fraction of the CMB rad i ation and emit black- 



body radiation at a higher temperature (iBabich et al 
intensity of the CMB radiation is then 

h{n) = [1 - r{fi)]B,[TcMB{n)] + T{n)B4Toort{n)], 



20071 ). The observed 



(3) 



where r is the optical depth and Tpprt is te mperature of the Oort Cloud objects 
along the direction n. iBabich et al.l (120071 ) constrained the mean, or monopole, 
contribution to the CMB spectrum using the Cosmic Backgr ound Explorer 



Far Infrared Absolu t e Spe ctrophotometer (COBE FIRAS) data (IFixsen et al. 



19971 : iMather et al.l . Il999l ) and placed upper limits on the total mass in the 
inner Oort Cloud for a variety of models of the size distribution of Oort Cloud 
objects. In this paper we will focus on effects that can make the Oort Cloud 
non-spherical and subsequently produce anisotropic spectral distortions in the 
CMB. 



We should emphasize that the induced spectral distortions considered in this 
paper are not the chemical potential [fi) or Compton-y spectral distortion that 
are commonly discussed in the literature. The signal produced in the outer So- 
lar System is that of a sum of weighted blackbodies at different temperatures, 
which in general does not have a blackbody frequency spectrum. The excep- 
tion is the Rayleigh- Jean's portion of the spectrum, where the weighted sum 
of blackbodies does produce a frequency spectrum described by a Rayleigh- 
Jean's distribution at a weighted temperature. This implies that low frequency 
data will not be very useful in allowing us to detect a signal originating from 
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the outer Solar System. See iBabich et al.l (120071 ) for a discussion of techniques 
to distinguish the spectral distortions produced in the outer Solar System from 
those caused by the CMB temperature anisotropics. 



The Oort Cloud is believed to be nearly spheri cal based on the statistics o f 
the inclination angles of new long period comets (IMarsden fc Sekaninal . Il97ll ). 
The observed apsides of these comets are not uniformly distributed across the 
sky, however the observed clusteri ng trends are believed to be due to the 



structure of the Galactic potential (jPelsemmd . 119871 : iDones et al.l . |2004| ). The 



Oort Cloud has been populated by objects that formed in the Sun's proto- 
planetary disc and were scattered multiple times by the four giant planets 
into increasingly higher eccentricity and higher energy orbits. This process 
continued until some external perturbation lifted the comets' perihelia well 
beyond Neptune's orbit, and then the comets became decoupled from the inner 
Solar System. The radial distribution of the Oort Cloud comets, which is set 
by how far the gas giant planets can scatter the comets before their perihelia is 
lifted safely outside the inner Solar System, strongly depends on the external 
environment of the young Solar System. The evolutionary processes that lifted 
the perihelia of comets beyond the planetary region of the Solar System must 
have also made shperical the distribution of comets which where initially in 
the Ecliptic Plane. Simulations of this formation process have found an inner 
edge of the Oort Cloud, defined as the semi-major axis at which the co mets 



are isotropica lly distributed, ne ar a = 3000 AU (IDuncan et al.l . 119871 ) and 
a = 5000 AU (IDones et al.l . 120071 ). It should be emphasized that objects with 



smaller semi- major axes exist, but may be more confined towards the Ecliptic 
Plane. This is the Scattered Disc which smoothly extends from the flattened 
Kuiper Belt to the spherical Oort Cloud. We will quantify the effect of the 
Scattered Disc by calculating the quadrupole moments of its thermal emission. 
This will exist irregardless of any externally induced asphericity in the Oort 
Cloud. 



An increasing body of evidence now suggests that the Solar System did not 

form in isolation, but in a stellar cluster which subsequently evaporated (j Adams fc Laughlin . 

200ll : iHester et al.l . 120041 ). This is not surprising since most low mass stars are 

observed to form in clusters that subsequently evaporate when the interstellar 

gas, which is required to gravitational bind the cluster , is expelled from th e 

cluster by either stellar winds or supernova feedback (ILada fc Ladal . l2003l ). 

Beyond this general trend, anomalous abundance ratios are observed in the 

Solar System that would indicate that a Type II supern ova occur somewher e 



in the vicinity of the Solar System during its infancy (iLooney et al.l . |2006| ) 



With a higher number density of stars in the cluster the rate of external per- 
turbations would have been higher and the spherical inner edge of the Oort 
Cloud would have extended to smaller semi-major axes. 



In this paper we will consider the dynamical effect that a star passing within 
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< 1 pc of the sun would have had on the Oort Cloud. And we calculate the sig- 
natures that would be produced in the CMB frequency spectrum. The outline 
of this paper is as follows. In §2] we describe how a non-sherical distribution 
of Oort Cloud objects may affect the measured brightness distribution of the 
CMB on the sky. In §2.11 we calculate the variance of the Oort Cloud signal 
due to Poisson fluctuations, whereas in §2.21 we describe the initial distribu- 
tion of Oort Cloud objects. In §2.31 we outline of method of calculation the 
perturbations in the orbital elements of the Oort Cloud objects, and in §2.41 
we calculate the distribution of impact parameters of stellar perturbers. In 
^ we derive an formula based on the coUisionless Boltzmann equation that 
allows us to analytically estimate the resulting quadrupole distortion from a 
perturbation. Our results are discussed in ^ In ^ we develop a Bayesian 
estimator that can be used to constrain the parameters of a stellar encounter 
from the observed CMB data. Finally §S] summarizes our main conclusions. 
Appendix |X] derives relationships between various measures the quadrupole 
moment and Appendix [B] describes the model of instrumental noise for the 
Planck satellite. 



2 Calculation 



We start by calculating how the non-sphericity of the Oort Cloud affects the 
measured angular distribution of the CMB intensity. The fluctuation in the 
CMB frequency spectrum along a given direction 

5h{h) = h{h)-B,{fcMB), (4) 

can be expressed as 

6Un) =Ar(n)^^ (5) 

+ /a f dxip{xi)—J—r[B^{TooTtixi)) - 5^(Tcmb)], 
J D[xiY 

where is the probability distribution of orbital elements describing the 

Keplerian orbits of the Oort Cloud objects around the Sun. The first term on 
the right is due to the extragalactic CMB temperature anisotropies. The CMB 
signal that comes from the surface of last scattering [z = 1000) has a black- 
body frequency distribution th at is parameteri zed by slightly different tem- 



peratures in different directions (jPodelsonl . l2003l ) . Once we have subtracted off 



the frequency spectrum corresponding to the average CMB temperature, Eq. 
(jl]), the CMB temperature anisotropies induce some spectral distortion since 
we have removed a blackbody at an incorrect temperature for that particular 
pixel. The second term is caused by the Oort Cloud objects in the outer Solar 
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System as described in ^ where D is the distance of the Oort Cloud object 
from the Sun. The absolute covering fraction is 



where R is the radius of the Oort Cloud object. The absolute covering fraction 
is related to the optical depth as r = /^/Z)^. Since the individual objects are 
optically thick, but sparsely fill the telescope beam function with which the 
CMB sky has been observed, the optical depth is simply related to the covering 
fraction of Oort Cloud objects. We add the term 'absolute' because fA is the 
covering fraction if all Oort Cloud objects were placed at a standard distance. 
We have implicitly assumed that the physical size and the orbital elements 
of the Oort Cloud objects are uncorrected. This assumption allows us to 
decouple the probability distribution of the orbital elements used in Eq. ([5]) 
from the size distribution, defined in Eq. ([7]) and used to calculate the absolute 
covering fraction. This is a reasonable assumption since the comets behave as 
test particles with respect to Sun and the perturbing stai[ll- For the internal 
density of the Oort Cloud objects we adopt a value of p = lgcm~^. 

In order to relate the mean optical depth to a total mass of the inner Oort 
Cloud, the differential mass function must be specified. We model it as a 
broken power- 

uT) 

(M) =AM-, M^i„ < M < Mbr, (7) 



dM 

The mass distribution yields the total mass through the integral 

dn 



/U/Tl 
dMM—{M). 



The empirical values of the above power-law indices are not known. For the 
Kuiper Belt, the slope of the high mass distribution has been restricted to the 



^ Note that it is not obvious that the initial distribution of orbital elements 
should be in dependent of the object's mass becau se non-gravitational forces, such 



as gas drag (IWeidenschiningj . 119771 : iRafikovl . 12004 ) and Pointying-Robertson drag 
(iBurns et al.l . Il979l ^. may have been important when the gas giant planets were 



scattering the comets into highly eccentric orbits. Also, coUisional effects, which 
are more important for sniall objects, may alter the format ion of the Oort Cloud 
(IStern k Weissmanl . [ioOll : ICharnoz Morbidenil . l2003l . [20oi ). These are issues that 



go beyond the scope of this paper. 

^ Typically the differential mass function is expressed in terms of comet radius, not 
mass. The power-law index with respect to radius Oradius is related to the values 
used in this paper as amass = (aradius + 2)/3. 
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range (3 = 2 ± 0.2 (jBernstein et all |2004J ). There are no strong observational 
constraints on the slope of the low mass distribution, a; ho wever, it cari be ca l- 
culated through the assumption of collisional equilibrium (jPan fc Saril . |2005| ) . 
The exact value of the slope then depends on the assumed material strength 
of the Kuiper Belt objects. Similiar constraints on the inner Oort Cloud mass 
distribution does not exist. Since collisions have not been numerous enough 
to affect the slope of the high mass distribution, it is reasonable to adopt the 
observed value of the slope of planetesimals in the proto-planetary disc and 
use it for the distribution of inner Oort Cloud objects. We will adopt = 2.16 
for the calculations in this paper. The collisional grinding of the Kuiper Belt 
objects, which justifies the above assumption of collisional equilibrium, is be- 
lieved to have occurred amongst the objects at their current location. The 
collision timescales in the Oort Cloud are much longer, so a state of collisional 
equilibrium cannot have been reach if the objects were always at their present 
locations. Fortunately collisions are believed to have been frequent during the 
formation phase of the Oort Cloud while the objects were still in the planetary 
region of t he Solar System (jStern fc Weissmaru . 1200 ll : ICharnoz fc Morbidellil . 
20031 . 12OO6I ). Again we will adopt the Kuiper Belt value of a = 1.83. Our choice 
of this value is conservative, since the primordial slope was likely steeper, pro- 
viding even better constraints than we will derive. 



We assume Mmax = ^Sedna = 6 X 10^^ g (jBrowru . 120071 ). as the upper mass 
limit. Our results are insensitive to changes in the upper mass limit for (3 > 2. 
We will take the minimum mass to correspond to objects of radius 1 mm. This 
quantity depends on the detailed formation history of the Oort Cloud and is, 
therefore difficult to calculate from first principles. The minimum object mass 
is partially degenerate with the total mass of the outer Solar System which is 
also unknown, so for the purposes of this paper we will fix the minimum radius 
at 1 mm. However, we note that the efficiency of blackbody emission is sup- 
pressed when the radiation wavelength is larger than the size of object emitting 
the radiation due to Kirchoff's Law (IGreenbergl . 1 19781 : iBackman fc Parescd . 
I993I : iLandau fc Lifshita . Il999l ). It should be emphasized that we only need 
our adopted model of the size distribution to convert a total mass of the 
Oort Cloud into an optical depth; if we simply choose a value for the optical 
depth then the rest of this paper would be independent of our assumed size 
distribution. 



For the location of the break in our mass model, we assume Mbr = 4 x 10^^ g 
(corresponding to i?br = 1 km). This value is established by requiring that 
objects of this size were likely to have undergone a collision during the time 
period of Oort Cloud's formation. The formation of the Oort Cloud most likely 
was a continuous process so it is difficult to precisely define its formation time. 
In simulations this time period can roughl y be estimate d asAtoort = 200 Myr 
based on a visual inspection (see Fig. 8 of iDones et al.l (120041 )). The collision 
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rate of objects of size i? in a protoplanetary disc is (IGoldreich et all |2004J ) 



(9) 



where a is the surface mass density of planetesimals, Q is the Keplerian 
velocity and p is the planetesimal density. We adopt minimum mass solar 
nebula value for the su rface density of condensates a = (20 AU/a)^/^ g cm~^ 
( IGoldreich et al.l . |200J) . Now the break radius can be calculated as 



R 



br 



crfiAtoort 

p 

1 km 

(a/20AU)= 



(10) 
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A full simulation would be required to precisely determine this quantity as well 
as its radial depen dence within the Port C loud. This is roughly consistent with 
the HST results of iBernstein et al.l (120041 ) for the Kuiper Belt. Our results are 
insensitive to the precise choice of these values since they are dominated by 
the smallest objects which have the largest surface area to volume ratio. 

The total mass in the Oort Cloud is calculate d by requir i ng the re to be IOMq, 
in objects _R > 1 km in the outer Oort Cloud ( Weissman . 1990l ). This figure is 
related to the observed number of new long period comets and implies there 
is approximately 10^^ comets present in the outer Oort Cloud. We use our 
adopted models to relate this number of comets to the total mass in the inner 
Oort Cloud. For simplicity we assume equal amounts of mass in the inner 
and outer Oort Clouds. Simulations of the formation of the Oort Cloud have 



implied that the inner Oort Cloud is more massive (IDuncan et al.Ul987h. while 



other s have implied that the outer Oort Cloud is more massive (iDones et al. 
20071 ). This uncertainity, as well as the uncertainity in other aspects of our 



model, should be remembered when interpreting our results. In Table [T] we 
present slightly different models and calculate the total number of objects, 
the total mass and the effective optical depths of these models. We quote the 
absolute covering fraction (/^i) divided by the mean squared distance 



/ daPia) — . 
J 



(12) 



The distribution of semi-major axes and the limits of integration are described 
in §2.2[ For definitiveness we adopt -P(a) oc a~^; Omax = 10^ AU and amin = 
10^ AU for our calculations in Table [H These parameter choices lead to the 
value 

— r- = 5.05 X 10"^AU"2_ (-^3) 

Next we calculate the microwave anisotropics induced by the thermal emission 
of the inner Oort Cloud. We consider the real-space quadrupole moments 
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i?min(cm) 


-Rbr(km) 


a 


^Oort 


Mtotai(Me) 


r = fA/{a^) 


0.1 


10 


1.83 


6.2 X 10^8 


1184 


1.8 X lO-'^ 


0.1 


1 


1.83 


1.0 X 10^9 


590 


2.9 X 10-^ 


0.1 


10 


2.0 


7.1 X 10^^ 


2565 


1.2 X 10~3 


0.01 


10 


1.83 


8.1 X 10^1 


1194 


2.1 X 10-5 



Table 1 

Total number (A'^oort); total mass (Mtotai) and absolute covering fraction (/a) of 
Oort Cloud objects for the adopted model in this paper (the first row) and slight 
changes to the model. Also shown is the optical depth (r) for each model. We 
calculate these quantities by requiring that 50M^ of comets {R> 1 km) exist. 



instead of the more conventional spherical harmonic moments, because it is 
easier to interpret the results in real space. These moments are defined as 



2 , . - - 1 . ^ 6I,An) 



Qij{v) = d h {hihj - -Sij)— , (14) 



CMS 

where Iqmb = -^^(Tcmb) is the homogeneous intensity of the CME0. We also 
quote our results as temperature fluctuations, which are more commonly used, 
in gi 

There are five independent quadrupole moments because the quadrupole ten- 
sor is symmetric and traceless. The Oort Cloud quadrupole moments are also 
frequency dependent. The standard CMB moments are frequency indepen- 
dentH because the usual sources of the cosmological temperature anisotropy, 
gravitational redshift and Thomson scatteri ng, do not depend on the frequency 



of the scattered radiation (jPodelsonl . |2003| ). Of course, we can relate the real 



space quadrupole moments to the spherical harmonic moments. These rela- 
tionships and the conversion between different measures of the quadrupole are 
discussed in Appendix \^ 



^ We should note that if the Oort Cloud was intrinsically spherical in a heliocen- 
tric coordinate system and the directional vectors used to measure the quadrupole 
moments, n in Eq. (jl4p . were defined in a geocentric coordinate system, then a 
small quadrupole moment would be observed. This quadrupole moments is roughly 
Q ~ {^(B/rooTt)'^^! ~ 10^^(5/ and therefore quite small. 

^ This statement is only true when the primordial CMB moments are expressed 
in terms on temperature, not intensity, flucuations. The Oort Cloud moments ex- 
pressed as either temperature or intensity fluctuations are frequency dependent. 
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2.1 Poisson Fluctuations 



The flux emitted by a single unresolved Oort Cloud object scales as 




At higher frequencies in the Wien portion of the spectrum the scaling has a 
steeper dependence on distance D. Consequently, even when we calculate the 
time average of the quadrupole moments 



the integrand still scales as an inverse power of i^. In the above equation we 
have used Kepler's second la-«0]to change integration variables from time (t) to 
true anomaly (/). In Eq. (fT6|) P is the orbital period of the Oort Cloud object. 
Therefore, an individual Oort Cloud object will produce the largest amount 
of signal per unit time when it is closest to perihelion and equivalently the 
comets that happen to be nearest to their perihelia will contribute the greatest 
to the overall signal per unit surface area. Since the comet's velocity is largest 
close to perihelion there will be fluctuations in the direction and amplitude 
of the quadrupole. These fluctuations will become larger at high frequencies 
where the signal's scaling with inverse distance is even stronger. The signal in 
a given direction is the average signal of all Oort Cloud objects weighted by 
their apparent angular size. If the relatively few nearby objects that appear 
large, have a higher temperature and move quickly, dominate the signal over 
the more numerous farther away objects, then the Poisson fluctuations will be 
larger. 

We adopted a continuum description of the Oort Cloud object masses and 
orbital elements, which implicitly assumed an infinite number of objects, when 
calculating the induced CMB spectral distortions in ^ In reality there are 
only a finite number of objects along any given line-of-sight and the signal will 
be dominated by the smaller number of objects that are closest to perihelion, 
so Poisson fluctuations will produce statistical fluctuations. These fluctuations 
will be estimated in this section. 

The surface brightness in a given direction, which is now a random variable, 
can be expressed as a sum over the contribution from any object which lies 

^ This states that a vector joining two masses will sweep out equal areas in equal 
times, which is equivalent to the conservation of the specific angular momentum, 




(16) 



(17) 



h = r^f. 
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along the line-of- sight, 



1 ^ 

6h{fi) = \2e'^^\R,/r,-\n-n,\) 
i=i 

X [i?,(T^,,J)-i?.(f,^b)], 



here hi is the position on the sky of the Oort Cloud object's centroid and 
the Ri is the Oort Cloud object's radius and its radial distance. The 6- 
function simply requires that the direction of interest intersects the relevant 
Oort Cloud object. Here N is the total number of objects in the sample that 
we are considering. We can directly calculate the Poisson fluctuations 

1 ^ 

{6h{n,)6h{n2)) = Y.i^^'HRi/D^ - - h,\) (19) 



X 



The Oort Cloud objects subtend a finite angle on the sky and Poisson fluctua- 
tions will cause CMB fluctuations if fii and n2 intersect the same object. The 
expectation value is taken over the assumed distribution of orbital elements 
and masses. 

We are interested in the intensity fluctuation averaged over some angular 
aperature, which we take to be a top-hat of angular size 9, 

-^^Un). (20) 

The power induced in the CMB by the Poisson fluctuations can be related the 
variance of the aperature averaged intensity fluctuations as 

'2^^^[BAf~^))^' ^^^^ 
where smoothing was done on the angular scale 9 = 27r/£. 

At a fixed time the average over the initial time of perihelon is equivalent to 
performing a time average over the orbital period which can be related to an 
average over the true anomaly. After we perform these averages we find that 
the power spectrum produced by Poisson fluctuations is 

d =— ^ / dMdadedeP{e)P{a)P{e)P{M) (22) 
27r A* J 



1 



cmb J 



The distributions of orbital elements and mass are given in Eqs. [23] & [TJ 
respectively. The integrals over the three orientation angles are unity. Although 
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/^codc 




oU LrJlZ 


1.6 X iU 


-32 


1.6 X iU 


-14 


100 GHz 


5.2 X 10" 


-32 


5.2 X 10^ 


-14 


217 GHz 


7.4 X 10" 


-31 


7.4 X 10" 


-13 


353 GHz 


2.5 X 10" 


-29 


2.5 X 10" 


-11 


545 GHz 


5.4 X 10" 


-27 


5.4 X 10" 


~9 



Table 2 

Fractional Poisson fluctuations in the number of Oort Cloud objects at different fre- 
quencies. The fractional fluctuation relevant for the calculations of our code (C™'^'^) 
and the real Oort Cloud (C^*^) are shown. 



we will apply this formula to the quadrupole = 2), it is applicable to the 
power spectrum on any angular scale. 

There are two distinct reasons that we are interested in Poisson fluctuations. 
The first involves understanding the accuracy of our numerical calculations, 
described in §2.31 We will perform our calculation with 10^ realizations of 10^ 
test particles, but due to the stronger weighting of the test particles closest to 
perihelia, the effective number of test particles is reduced and the fractional 
fluctuation amplitude increased. The second is related to the statistical signif- 
icance of any observed CMB spectral distortion. As Oort Cloud objects move 
along their orbits the closest objects, and therefore the objects that produce 
the largest signal, will change. These Poisson fluctuations could potentially 
limit the statistical signiflcance and therefore usefulness of any measurement 
of CMB spectral distortions. The results of our calculations are shown in Table 
[2] for the fractional fluctuations relevant to both our code and the real Oort 
Cloud. 

We conclude that statistical fluctuations will not be important when interpret- 
ing the observational data if our model of the size distribution of Oort Cloud 
objects is reasonable. These fluctuations will be much smaller than standard 
experimental noise. However, the statistical fluctuations present in our numer- 
ical calculations, which contain signiflcantly fewer test particles than the true 
number of Oort Cloud objects, will be considerably larger and comparable to 
the magnitude of the effect that we are interested in. Therefore we will cal- 
culate the distribution of quadrupole moments by calculating the quadrupole 
at 10'^ different times. From this distribution we can estimate the magnitude 
of the perturbation on the shape of the Oort Cloud, as well as, the effect of 
Poisson fluctuations on the resultant CMB spectral distortions. 
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2.2 Initial Conditions 



The comets in the Oort Cloud, in the absense of external perturbations, ef- 
fectively act as test particles that move on Keplerian orbits around the Sun. 
The hypothetical stellar encounter simply perturbs the comets onto Keplerian 
orbits with new orbital element J^. The orbits are completely described by six 
orbital elements - (i) the semi-major axis (a); (ii) the eccentricity (e); (iii) 
the inclination of the orbital plane (iv) the longitude of perihelion (cj); 
(v) the longitude of the ascending node and (vi) the time of periheliorLLl 
(0). The inclination angle, the longitude of perihelion and the longitude of the 
ascending node, are rotation angles used to convert the coordinate system of 
invariant plane (which is different for every comet) to a standard reference co- 
ordinate system. With these orbital elements the motion of the comet around 
the Sun is completely specified. 

We set the initial conditions of the Oort Cloud objects by randomly sampling 
the uniform distributions 

50 AU 

0.3 < e < 1 , (23) 

a 

-l<cos/<l, (24) 
< < 27r, (25) 
< < 27r, (26) 
< e < P, (27) 

where P = / [GMqY^"^ is the test particle's orbital period. The maxi- 

mum eccentricity is determined by the condition that the comet's perihelion 
lies far enough away from Neptune that it becomes dynamically decoupled. 
This condition roughly requires Tp = a{l — e) > 50 AU, so the maximum ec- 
centricity depends on on the semi-major axis of the orbit. Simulations of the 
formation of the Oort Cloud have found that the probability distribution of 
semi-major axes has a power-law form P{a) oc g ""^, but have d i sagre ed on 



the value of its index 7. Originally it was stated by lDuncan et al.l (119871 ) that 



7 ~ 3.5, but updated simulations with more re alistic initial c ondit ions for the 



comets have found a shallower value, 7 ~ 1.5 (IDones et al.l . 120071 ). When we 
present numerical results in §l]we adopt the value of 7 = 3. We will assume 
that the inner edge of the Oort Cloud is a = 10'^ AU. We take the inner edge 
to be semi-major axis beyond which the Oort Cloud becomes roughly spheri- 
cal. The Scattered Disc lies inside of this radius. And we pick the outer edge 
to lie at a = 10^ AU. We then randomly choose an initial time < ti < P 



^ We also do allow the comets to be ejected from the Solar System if the final 
eccentricity if greater than one. 

^ Ordinarily r is used to designate the time of perihelion passage, we chose a dif- 
ferent symbol to prevent confusion with the use of r for the optical depth. 
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to calculate the initial position in order to evaluate the above perturbation 
equations. Then we choose another random final time < tf < P in order to 
evaluate the final positions of the comets. 

The initial conditions of the Oort Cloud before the stellar encounter depend 
on the rate at which the comets are scattered to high inclination angleo. At 
smaller distances, the isotropization processes are less efficient and the cloud 
will be smoothly compressed into a torus that extends out of the ecliptic plane. 
The initial conditions depend on the efficiency of these relaxation processes. 
These processes include the effects from passing stars and molecular clouds, the 
tides of the Galactic disc, and the potential of the parent molecular cloud inside 
of which the Sun had formed. Placing a large amount of mass into the non- 
spherical Scattered Disc or any primordial asphericity in the Oort cloud will 
only increase the anisotropic spectral distortions. The intensity fiuctuations 
induced by the Scattered Disc estimated in this paper can be viewed as a lower 
limit for the signal. 

We estimate of the infiuence of the initial conditions by calculating the ex- 
pected quadrupole signal from the Scattered Disc. We assume that the total 
mass is O.IM^ and the semi-major axes lie between 10^ AU and 10'^ AU. Most 
importantly we choose the disc's opening angle to be 30 degrees. The six 
quadrupole moments at three frequencies, u = 30 GHz (black), u = 143 GHz 
(red) and z/ = 217 GHz (blue), are shown in Fig.[Tl The mixed quadrupole mo- 
ments, such as Qxy, all vanish within Poisson fiuctuations. As expected Pois- 
son fiuctuations become larger at higher frequencies. The diagonal quadrupole 
moments, such as Qxx, do not vanish. The Q^z moment does not vanish be- 
cause the disc is compressed into the Ecliptic Plane. The traceless condition, 
Qxx + Qyy + Qzz = 0, forccs the other two diagonal moments to be non-zero. 
Since the scattered disc is azimuthally symmetric by construction, Q^x and 
Qyy are forced to be equal. 

2.3 Orbital Perturbation 

The external perturbation from a passing star acts both on the Oort cloud 
as well as the Sun which are bound together, and so only the relative (tidal) 
acceleration between the comet and the Sun is relevant 

Aa =aoort — o-Q, 

GM,(h - r) GM,h 
- |b-r|3 6^' (2^) 



^ The initial distribution of oj and should have been roughly uniform even in the 
proto-planetary disc. 
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Fig. 1. Histograms of the six quadrupole moments at three different frequencies 
u = 30 GHz (black), u = 143 GHz (red) and u = 217 GHz (bkie). The mixed 
quadrupole moments all vanish within Poisson fluctuations and the diagonal mo- 
ments have are nonzero at a statistically significant level. 

where Mp is the mass of the perturber and b is the vector from the Sun to the 
perturber. 

Kepler's third law, T oc a^/^, implies that an object on an orbit with a = 
1000 AU has an orbital period, T = 10^^ s. If the impact parameter of the per- 
turber is 6 ~ 2000 AU = 0.01 pc and it is traveling at a speed of f = 30 km/ s, 
then the effective duration of the impulse would be Tpert ^ 2b/ v = 2 x 10^° s. 
The orbital period is typically longer than the duration of the impulse, and so 
we will use the impulse approximation for simplicity. Numerical experiments 
have found that the validity o f the impulse approxima tion even extends to the 
regime where T 0.2 x Tpert (lAguilar fc Whitd . Il985l ). For test particles with 
larger semi-major axes the accuracy of the impulse approximation becomes 
increasingly better. 



Clearly, ge through the Oort Cloud would strongly pe r turb t he comets 

that lie nearly along the perturber's trajectory (jWeissmanl . Il990l ). Such an 
effect is believed to be responsible for the comet showers that occasionally take 
place (as the strongly perturbed comets either enter the inner Solar System 
or get ejected from the Oort Cloud), but will these strongly perturbed Oort 
Cloud objects will have little effect on results because they are assumed to have 
been removed from the Oort Cloud. In our calculation will assume that all test 
particles that are more strongly attracted to the perturber than the Sun are 
ejected. Initially this has a large effect on the angular structure of the Oort 
Cloud, however the effect is quickly diminished due to phase mixing. After a 
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test particle is ejected from the Solar System a hole is left in the Oort Cloud 
that orbits around the sun with the orbital parameters of the original test 
particle. This results in an underdensity of test particles that would produce 
anisotropy in the same manner as an overdensity would. However the large 
spread in orbital periods, due to the large spread in semi-major axes, means 
that the signal vanishes over a timescale corresponding to the orbital period 
of the test particles. Of course, the mean signal will remain the same as this 
only depends on the total covering fraction of Oort Cloud objects. However 
the quadrupole, or any angular moments, will vanish as the holes left by the 
ejected test particles lose spatial coherence. 



In the impulse approximation, only the component of the force p erpendicular 



to th e perturber's motion remains after its path is integrated over (IBinney fc Tremaind . 



19871 ). Therefore we must apply the projection operator 



vv, (29) 



to the tidal force, where v is the velocity of the perturber. The components of 
the projected tidal force, written in a coordinate basis attached to the comet, 
are 



R = Br-T ■ Aa, (30) 
T = eTTAa, (31) 
N = Bn ■ T ■ Aa, (32) 



here e.^ is the radial basis vector, is the tangential basis vector and is the 
normal basis vector. The Lagrange equations describe how the perturbation 
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changes the orbital elements ( Danby . 1962 ; Murray fc Dermott , 2000l ) 

da 2a3/2 



dt 

de 
'dt 
di 
dt 

dn 
Itt 

duj 
Itt 



dB 
'dt 



Ml-e^ 

/a(l-e2 
/i 



:[R esin/ + T(l + esin/)], 



D • ^ , T-r ^ , e + cos / 

sm / + T(cos / + — -) 

1 + e cos J 



/a(l-e2) A^cos {00 + f) 
1 + e cos / 



/a(l-e2) A^sin(cj + /) 
/i sin?(l + ecos/) 



(33) 

(34) 
(35) 
(36) 



la 1 



■[— -Rcos/ + Tsin/ 



2 + e cos / 
1 + e cos / 



■ COS? 



dt 

m{Q - 1) 



Ml 



-e sin / 



(37) 
(38) 



+ R 



a^(l-e^) 



+ 3T(e-t), 



/i(l 



2 cos / , 

1 + e cos / e ' 

n 

—e sin / 



T 



a2(l-e2)sin/(2 + ecos/) 



e/i 



1 + e cos / 



where /i = GMq. 



Within the impulse approximation, we assume that the internal configuration 
of the solar system is frozen during the close encounter with the perturber; it 
is this assumption that causes the force parallel to the perturber's velocity to 
exact cancel. Therefore the change in the orbital element, x, can be expressed 
as 



dx2h 
Ax ^ — . 

dt V 



(39) 



The new distribution of orbital elements will be solely due to the addition of 
velocity perturbations, because all previous comets with the given elements 
will have been perturbed to other orbits 



Pi{xi}) = J Po{xi - /\Xi)D{Axi) 



(40) 



In the test particle limit, the mass of the comet should not affect how it is 
perturbed by the stellar companion. Therefore the distribution function of 
the Oort Could comets should be factorizable into the distribution of comet 
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orbital elements and the mass function of comets 



foon = j^{M)xp{{x,}). (41) 



We will use the orbital element distribution function in Eq. (HTj) to calculate 
the CMB intensity fluctuation using Eq. In principle the orbital element 
distribution function can be used to construct a distribution function param- 
eterized by the Keplerian integrals of motion. Since the distribution function 
is six dimensional and we are interested in observable signatures in the CMB, 
we will not focus on this aspect of the problem. However we will mention a 
condition that the distribution function must satisfied if an angular signal will 
be produced. For a quadrupole, or any angular moment, intensity fluctuation 
to be non-zero, the distribution of Euler angles {Q,u,i) must be correlated 
with the distribution of semi-major axes and eccentricities. In our calcula- 
tions we initially assume that all the orbital elements are uncorrelated and 
that the Oort Cloud is statistically isotropic. The stellar perturbation induces 
the necessary correlations to produce a statistically significant signal. 



2.4 Likely Encounter Parameters 



In addition to determining the initial conditions of the Oort Cloud objects we 
need to specify the properties of the perturber. Observations of local stars can 
identify their impac t parameter with respect to the S olar System only up to 



10 Myr in the past (IGarcfa-Sanchez et al.l . Il999l . I2OOII ). Of course, we can use 



the initial mass function of stars and their velocity dispersion with respect 
to the local standard of rest to make statistical statements about the likely 
properties of the perturber, however for a single perturber that gets to the 
closest approach and has the biggest effect, such statistical statements are not 
very useful. In ^ we will discuss a Bayesian method that allow us to use a 
detected signal to statistically constrain the properties of the perturber. Given 
a number density and a relative velocity with respect to the Sun the minimum 
impact parameter of the perturber, Rmin, can be estimated by requiring that 
the probability of encounter during some time period. At, is of order unity 

P = 7TR^n,t^,,vAt ^ 1. (42) 



The environment of the Solar System might have been quite different in the 
past than it is today. As mentioned in ^ the Sun might have formed in 
a stellar cluster and so the number density of neighboring stars may have 
been much higher. Also as the Solar System undergoes epicyclic and vertical 
oscillations with respect to the Galactic disc or passes through spiral arms, the 
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Model n(pc-3) Vreii'km s'^) At(Gyr) i2^i„(AU) 



Present Neighborhood - 


A 


0.1 


30 


4 


1100 


Present Neighborhood - 


B 


0.1 


30 


1 


2100 


Present Neighborhood - 


C 


0.1 


30 


0.5 


3000 


Cluster - A 




1000 


1 


0.01 


1100 


Cluster - B 




100 


1 


0.01 


3700 



Table 3 

Characteristic number densities, velocities, time intervals and minimum impact 
parameters for several models of the local Solar neighborhood. We adopt the pa- 
rameters of 'Present Neighborhood - B' for our calculations. 



environment also changes. The closest encounter probabilities can be adjusted 
to account for the changing Solar environment. 



A variety of relaxation processes - pa. s sing s tars and giant molecular clouds 
( jPuncan et al.l . 119871 : iHut fc Tremaind. Il985[). pr ecession due to the potential 
of the Galac tic disc (jHeisler fc Tremaind . Il986l ) and the self gravity of the 



Oort Cloud (ITremaind . l2005l ) - can isotropize the orbits of objects in Oort 
Cloud and reduce the non-sphericity produced by the stellar encounter. Of 
course the possibility of eliminating these signatures does not imply that the 
Oort Cloud will be spherical today, but that it will possess the non-sphericity 
induced by the closest stellar encounter that occurred at a time less than 
the appropriate relaxation timescale. The self-gravity of the Oort Cloud was 
determined to pr oduce the largest effect; the calculated timescale was 1 Gyr 
( Tremaind . boosi ) . More calculations are needed to determine the effect on 
nearly parabolic orbits, which spend considerable amounts of time at different 
distances and therefore experience different gravitational perturbations. We 
will leave a detailed analysis of these relaxation processes to future work. 
Therefore we will consider the minimum impact parameters for the current 
Solar neighborhood with a variety of time intervals - At = 0.5, 1, 4 Gyr. These 
models are described in Table [31 



Of course rare events (P <^ 1) sometimes occur and the above figures should 
be seen as estimates of the mean of the distribution from which the actual 
closest impact parameter has been randomly drawn. The probability that the 
closest encounter occurs with impact parameter R is equivalent to the joint 
probability that one encounter occurred at R, P{1\R), and no encounters 
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Fig. 2. Plot of the probability distribution of minimum impact parameters for the 
models described in Table ([3]). The curves in a left-to-right manner correspond to 
Neighborhood- A (black, solid curve), Cluster- A (green, long dashed curve), Neigh- 
borhood-B (red, dotted curve), Neighborhood-C (blue, dashed curve) and Cluster-B 
(cyan, dot-dashed curve). We adopt the Neighborhood-B model for calculations in 
this paper. 

occurred with r < R, P{0\r < R), 

P{Rmin) =PiO\r<R) X P(l|i?), 

= /^p(r)rfre-""'™^* x {7iR'^nvAt)e-''^'''''^\ 
Jo 

g-TrR^niiAt ^ ^-wR'^nvAt'^ (43) 

This probability has a maximum at P = ln2 ~ 0.7, instead of P = 1 in Eq. 
(1421) . The probability distribution of Rmin for the variety of environments are 
shown in Fig. [2l 



3 Analytic Estimates 



Before we present numerical results in §1] it is worthwhile to look at order- 
of-magnitude estimates of the effect. Using the impulse approximation it is 
straightforward to determine how the perturbation affects the test particles. 
The velocity and kinetic energy are directly changed by the impulse, but the 
potential energy remains the same since the distance from the Sun has not 
changed. If the system was self-gravitating and the traditional virial theo- 
rem applicable, the system would relax to a new equilibrium distribution and 
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the moments of the gravitational p otential tensor would direc tly reflect the 
changed velocity dispersion tensor (IBinney fc Tremaind . Il987l ). Relating the 
velocity or energy kick to a quadrupole moment requires a generalization of 
the commonly used virial theorem as we show below. 



3. 1 Virial Equations 



The ensemble of test particles in the Oort Cloud satisfies the collisionless 
Boltzmann equation 

1^ + ■ V/ - V<l> • V./ = 0, (44) 

where the distribution function depends on position, velocity and time, al- 
though we will suppress this functional dependence for brevity. We will assume 
that the Oort Cloud is not self-gravitating, so the potential through which the 
test particles are moving is solely determined by the Sun's gravity 

$ a; = -. 45 

r 

Assuming the distribution is time- independent, which is equivalent to assum- 
ing that Poisson fiuctuations in the final moments are negligible (see §2.ip . we 
find the constraint equation 

GMo'^^ = v-Vf. (46) 

In simple words this equation states that the number of test particles ad- 
vected from a phase space volume element equals the number of test particles 
accelerated into this volume element by the Sun's gravity. If we were deriving 
the traditional virial theorem we would take velocity and position moments 
of this equation or the time-dependent version of it, Eq. (jS]). However we are 
interested in relating the velocity dispersion tensor to the intensity quadrupole 
moments, not the potential energy tensor. 

The quadrupole moment can be expressed as an integral over the Oort Cloud 
volume 

f Ia 1 

Qij{v)By{TcMB) = J d^rd^vf-^{xiXj - -r^6ij) 

X [5,(Toort)-5,(TcMB)]. (47) 
The optical depth and the distribution function are defined such that 

J d^rd^vf{x,v) = 1. (48) 
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The blackbody frequency spectrum is in general a complex function of Toort, 
although at low frequencies it scales as i?i,(Toort) oc Toort- In order to de- 
rive analytic results we will express the blackbody spectrum's dependence on 
distance as a power law, 

5.(Toort) = (49) 



where the exponent is defined as 



P 



1 din B^{T) 

2 d\nT ' 



(50) 



This exponent is only well defined in Rayleigh- Jean's limit where it equals a 
constant, f3 = 1/2. At higher frequencies its value grows rapidly and is strongly 
frequency dependent. Below we will use (3 = 1/2, but comment on how the 
higher frequency cases will affect our results. 



We first multiply Eq. (H6|) by vxr^^^^^\ Now following the basic steps in 
deriving the virial theorem we will integrate this equation over phase space 
and then integrate by parts to eliminate the position and velocity derivatives 
acting on the distribution function. This will allow us to relate the first part, 
the XiXj term, of the quadrupole moment to the velocity dispersion tensor. 
The second part, the —Sij/3 term, can be related to the velocity through the 
energy equation of a Keplerian orbit 



E 



GM, 







GMq 
2a 



(51) 



For simplicity, here and in the next subsection we assume zero eccentricity 
orbits with = GMq/q. 



The final equation relating the quadrupole moments to the velocity distribu- 
tion of the test particles is 



Q,,(z/)5,(Tcmb) = /a/ d?rd'vf{x,v){A, 
(2 + (3){xiVj + ViXj){x ■ v) 



ViVj 



GMqT^+P 



2GM0r3+/3 



VjVj 

GMr^r 



{XiVj + ViXj){x ■ v) 



GMqV^ 







3 


[gMq 




( 


3 


[gMq 



2+0 



(52) 



3. 2 Perturbations 



In this subsection we will make several approximations that will allow us to 
calculate analytic expressions for the quadrupole moments. We will adopt 
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the impulse approximation as well as the tidal approximation for the relative 
acceleration of the test particles with respect to the Sun 



Aarei = ^[r-3bir-b)]. (53) 

In this expression for the relative acceleration, Mp is the perturber mass, r 
is the radial vector between the Sun and the Oort Cloud object, and b is 
the impact parameter vector between the Sun and the perturber. We are able 
to use the tidal approximation in this section because we will assume a thin 
shell distribution of semi-major axes for all the test-particles. In our numerical 
calculations we assume a power-law distribution of semi-major axes between 
a = lO^AU and 10^ AU. With this latter distribution our analytic model 
assumption of the tidal acceleration is invalid because for the appropriate 
values of -Rmin the perturber passes through the Oort Cloud. 

With these assumptions for the analytic model, the velocity kick is, 

Av = -^^{I - vPvP)[r - 3b{r ■ b)], (54) 

where we project out the component of the force parallel to the perturber's 
trajectory. For later convenience we will define 

^ = n-3b{h-b) -vP{n-vP), (55) 

such that 

A?? = - 



2GM„r ^ 

Av = ^ ^. (56) 



Here we have assumed 

Vrel ■b=0, (57) 

which directly follows from our use of the impulse approximation. Working 
within the impulse approximation the direction of the impact parameter is 
almost entirely determined by the direction of the encounter velocity. This 
constraint from the impulse approximation allows us to simplify the form of 
the above velocity kick. There is still phase information that is needed to 
completely specify the geometry, Eq. (IFTI) is simply one equation for two un- 
knowns. For example, the direction of the encounter only determines a plane 
within which the direction of the encounter velocity must lie; the azimuthal 
angle of the encounter velocity must still be specified. Equivalently, the di- 
rection of the encounter velocity determines a great circle on the sky along 
which the encounter must have occured; the longitudinal position along that 
great circle must still be specified. The azimuthal angle of the relative velocity 
should have uniform distributions between and vr (recall the v — > —v de- 
generacy mentioned above), whereas the longitudal angle that determines the 
direction of the impact parameter should be uniformly distributed between 
and 271. 
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There is no change in the velocity dispersion, or the quadrupole moments, at 
first order in Av. This is a consequence of the fact that we assumed the Oort 
Cloud to be statistically time-independent and axisymmetric. At second order 
the velocity kick will produce a non-zero quadrupole moment. 



Q.A^)B^{Tcmb) - ^^^^ ^ 

, (2+/?) 



- B,{TcMB)[^i^j - i^iUj + ^jni){h ■ ^) - -6,j^ ■ *]}, (58) 

where we have assumed that the Oort Cloud is made of a thin shell of zero 
eccentricity test particles all with semi-major axis a. Since we assume that all 
of the Oort Cloud objects lie at a single semi-major axis a we can replace the 
absolute covering fraction defined in Eq. |6l with the optical depth r. 



The quadrupole tensor can now be directly evaluated as 

BuiTooTt) 



16nGM^ra^ 



156^2 



vel 



_B^(Tcmb) 

where the configuration tensor is defined as 



(59) 



AT -L ORB 

= [(3 - 0)v^v, + 3(1 + 3/3)6,6,] + —^^5,,. 



(60) 



Note that Vi and hi are components of b and v^^ which have unit norm. 

Taking some characteristic numbers, such as Mp = IM©, [3 = 1/20 
30 km s-\ a = 10^ AU and 6 = 2 x 10^ AU, we find 



Q{u) = 3.9 X 10" 



10-6 



g.(lOK) 
5,(3 K) 



(61) 



Ignoring the directional dependence described above, the amplitude of Q varies 
from 1.4 X 10"^ at low frequency to 2.3 x 10"^ at u = 300 GHz and 8.9 x 10"^ 
at z/ = 600 GHz. Our choice of r = 10"^ is close to the value of our adopted 
model in Table [H Since we have assumed that initial distributions of the test 
particle sizes and orbital elements are independent, we can simply scale the 
value of Q(z/) to the appropriate value of r. 



^ For simplicity we choose to evaluate the configuration tensor at (3 = 1/2 for all 
frequencies despite the fact that /3 and therefore Q will grow at higher frequencies. 
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Fig. 3. Histogram of the induced CMB quadrupoles at = 30 GHz before (red 
curve) and after (black curve) the perturbation. The quadrupole moments are in 
units of mean CMB intensity. 

4 Numerical Results 

Next we present our numerical results. We start by randomly choosing the ini- 
tial orbital elements of 10^ Oort Cloud objects according to the distributions in 
Eqs. (1231) . Then for a given trajectory and mass of the external perturber, we 
use Eqs. (1331) to map the initial orbital elements to the final orbital elements. 
The resulting quadrupole moments, calculated using Eq. f|T^ , are determined 
for 10^ distinct initial and final times in order to estimate the effect of statis- 
tical fluctuations on the outcome and therefore determine if the perturbation 
produced a statistically significant change in the structure of the Oort Cloud. 

These quadrupole distributions are shown in Fig. [3] for z/ = 30 GHz, Fig. H] 
for V = 353 GHz and Fig. [5] for u = 545 GHz. We choose these frequencies 
to correspond to the frequency bins of the Planck satellite. The statistical 
fluctuations are greater at high frequency because the induced signal scales 
as a stronger power of radial distance, as seen in Eq. f|T5l) . For this particular 
example the direction of impact parameter was B = (0.909, 0, —0.416) and 
the trajectory of the perturber was Vp = (0, 1,0). The perturber's mass was 
Mp = IMq, its relative velocity Vp = 30 km s~^ and impact parameter b = 
2100 AU. 

From the distributions of the quadrupoles in Figs. (jSHS]), we can determine that 
the non-zero value of the induced quadrupoles is statistically significant. Al- 
though our numerical calculations suffer from Poisson fluctuations, the signal 
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Fig. 4. Same as Fig. ([3]), but at = 353 GHz. 





Fig. 5. Same as Fig. ([3]), but at = 545 GHz. 

from the real Oort Cloud will not exhibit the relatively large changes over an 
orbital period if our model for the size distribution of Oort Cloud objects and 
estimate of the total mass are reasonable. The quadrupole moments shown in 
Figs. (l3]-[5l) can be converted to a value of C2, the Fourier space quadrupole 
power spectrum, using the formulae derived in Appendix \M The black points 
in Fig. [6] are the quadrupole power spectrum, C^=2, as a function of frequency 
along with the expected Planck instrumental noise errorbarg I. The expected 
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The instrumental error bars appear anomalously low because we are accustomed 
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Fig. 6. The quadrupole power, C^=2, at several different frequencies. The power 
spectrum is dimensionless since the Oort Cloud quadrupole moments are expressed 
in terms of the mean CMB intensity. The black points are the Oort Cloud power 
spectrum with error bars appropriate for Planck instrument noise. The blue points 
are the power spectrum produced by our model of the Scattered Disc. The red points 
are the theoretical power spectrum calculated with the best fit WMAP cosmology, 
corresponding to C2 = 1.76 x 10"^'^ in terms of temperature fluctuations, with 
cosmic variance error bars. Since we chose to represent the fluctuations in terms 
of intensity and not temperature, the primordial temperature anisotropies do not 
produce a frequency independent power spectrum. 

statistical fluctuations due to Poisson fluctuations in the number of observed 
Oort Cloud objects are much smaller for the models considered in this paper 
(see Table E]). Of course, the observability of the signal will primarily be af- 
fected by all of the astrophysical sources that emit non-blackbody radiation, 
such as the standard Galactic foregrounds, thermal Sunyaev-Zel'dovich emis- 
sion, etc. The signal originating in our Solar System must be observed against 
this background. The full calculation of the power spectrum as function of 
frequency of these various backgrounds is beyond the scope of this paper. 

Also shown in Fig. [6] is the quadrupole power spectrum of the Scattered Disc 
(blue points) and the theoretical power spectrum predicted by the WMAP 
cosmology model (red points). The theoretical power spectrum is frequency 
independent when expressed in terms of temperat ure fluctuations and has a 
value of C2 = 1.76 x 10~^° ( iHinshaw et al.l . |2006| ). Since we plot the inten- 
sity fluctuations, and not the temperature fluctuations, the theoretical power 
spectrum displays a slight frequency dependence. 



to seeing the cosmic variance errors bar dominate at low i. 
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Our assumption that the distributions of test particle mass and orbital ele- 
ments are independent allows us to vertically scale the curves in Fig. [6] as 
we change The predictions for the power spectrum scale as f\. Although 
we have separately plotted the quadrupole power spectra for the Oort Cloud 
and the Scattered Disc, in reality we will observe the combination of the two. 
We chose to separate the two contributions in Fig. [6] because the Scattered 
Disc is intrinsically aspherical and the Oort Cloud's asphericity is induced by 
a perturbation. Even if the Oort Cloud was perfectly spherical, the scattered 
disc, which is believed to smoothly compress from a spherical shape at the 
inner edge of the Oort Cloud into the Ecliptic Plane at the outer edge of the 
Kuiper Belt, will produce a minimum quadrupole moment. 

While we calculate the corresponding quadrupole power spectrum values for 
comparision with CMB values, we should emphasize that the individual quadrupole 
moments of the induced spectral distortions contain phase information. It is 
assumed that the CMB temperature anisotropies constitute a Gaussian ran- 
dom field and the observed a^^'s are interpreted as being a random realization 
drawn from a Gaussian probability distribution function with variance Ce. In 
which case, the only meaningful information in a set of CMB data are the mea- 
sured values of the Ces. For the anisotropic spectral distortions induced by the 
Oort Cloud the particular values of the individual a2m's, or equivalently the 
Qij (see Appendix |X] for relations that link these quantities), contain impor- 
tant information about the structure of the Oort Cloud. This is an important 
distinction to remember and one that we will explore in more detail in §3 

The instrumental noise error bars appear anomalously small in Fig. O Typi- 
cally on large scales cosmic variance dominates instrumental noise; for Planck 
the power spectrum should be cosmic variance limited until i ~ 1500. The 
error bars on the theoretical primordial quadrupole power spectrum in Fig. 
[6] are due to cosmic variance. Although the Oort Cloud and Scattered Disc 
power spectra are never distinguishable from the theoretical power spectrum 
at a single frequency, we can use the fact the primoridal CMB tempera- 
ture anisotropies are frequency independent to utilize the low frequency data 
(where the Oort Cloud signal is effectively zero) to eliminate the cosmic 
variance uncertainity in the primordial temperature anisotropies at high fre- 
quency. Namely the low frequency temperature anisotropies can be used to 
exactly predict the high frequency temperature anisotropieJ"!. The fact that 
cosmic variance limits our ability to measure the true power spectrum is irrele- 
vant because we only care about the particular realization of the temperature 
anisotropies in our Universe. This should allow us to use the temperature 
anisotropies in a single frequency bin to remove them from the other frequen- 
cies and detect the Oort Cloud signal whenever it is larger than instrumental 



This statement depends on the level of foreground contamination in the low fre- 
quency bands 
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noise. 



5 Bayesian Analysis 



The observed CMB can be used to constrain the parameters of the stellar 
encounter if a quadrupole, or any higher order moment, spectral distortion 
were to be detected. The posterior likelihood of the Oort Cloud and external 
perturbation parameters, collectively designated as T, given the observed data, 
designated as is 

C{T\D) (X C{D\T)C{T), (62) 

where we have used Bayes' Theorem to invert the conditional probabilities. 
The external perturbation parameters depend on the mass of the perturber, 
Mp] the impact parameter of the encounter, b; the velocity of the perturber 
with respect to the Solar System, Yrei] and the absolute covering fraction 
of Oort Cloud objects, There exist degeneracies between some of these 
parameters. For example, the projection matrix, defined in Eq. fl29|) . is invariant 
under the change v — > —v since the force parallel to the trajectory cancels 
in the impulse approximation, and the total impulse only depends on the 
combination Mp/b'^Vrei- Plots of the parameter constraints will directly display 
these degeneracies as the data is only sensitive to particular combinations, 
although the a priori probabilities will constrain these parameters to some 
extent. 



The likelihood function of the data given the stellar encounter parameter can 
be expressed as a functional integration over the true CMB signal (S), which 
will be a combination of primordial temperature anisotropics and Oort Cloud 



emission, 



£{D\T) = J dS/:{D\S)£{S\T), 



(63) 



where the likeliho od of t he data given the tru e CMB signal is simply set by 
instrument noise (IKnoxl . 1 19951 : iDodelsonl . |2003| ) 



C{D\S) 



(d-a)tN-i(d-a)/2 



V2tt det 



(64) 



where N is the noise covariance matrix, d is the observed signal and a is the 
true underlying CMB signal. The observed signal is not uniquely determined 
by the underlying theory which only specifies the probability distribution func- 
tion of which the signal is a random realization. The likelihood function of the 
CMB temperature anisotropics is 



C{S\T) 



1 

V2tt detC^ 



-atC-ia/2 



(65) 
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where C is the signal covariance matrix. The signal covariance matrix is de- 
fined as the expectation value of C = (aa^), with a similar definition for the 
noise covariance matrix N. The likelihood of the Oort Cloud signal given the 
theory is set a probability distribution is determined by Poisson fiuctuations 
as discussed above in ^2.11 



The a priori probability of Oort Cloud and stellar perturbation parameters 
should not be modeled as the product of the individual parameter distribu- 
tions. Observations of the local Solar neighborhood have found that different 
stellar spe ctral types have different number densities and different velocity 
ellipsoids ( iBinney fc Merrifieldl . 1 19981 : iKing et al.l . Il990l ) . We should be mind- 
ful that the Sun's local environment might have been quite different in the 
past as the Sun moves through different Galactic environments. The a priori 
distribution can be written in terms of conditional probabilities as 



CiD)=p{b\Vrel,M) pib\v) p(v,e/|M) p(M) 



(66) 



where p(M) is the distribution of stellar masses, averaged according to stellar 
evolutionary theory as 

p(M) = vl/(M)^, (67) 

where \1/(M) is the appropriate initial mass function, t{M) is the lifetime of 
stars of mass M and At is the lifetime of the Solar System. The distribution of 
veloc ities is assumed to follow a Schwartzschild distribution (IBinney fc Tremaind . 
1987h 



with t he velocity ellipsoid (S) dep ending on the stellar spectral type ( iKing et al.l . 
1990l : IBinney &: Merrifieldl . Il998l ). As discussed in ^ the direction of the im- 
pact parameter is the determined by the direction of the encounter velocity 
up to a phase angle. The magnitude of the impact parameter depends on the 
number density of stars of a given mass and their relative velocity with respect 
to the Su Given these parameters, the distribution of impact parameters 
can be calculated as outlined in §2.41 The a priori probability of J'a is com- 
pletely unknown, except that an upp er limit can b e set u sing constraint on 
the mean CMB spectrum distortion (IBabich et al.l . 120071 ). We can consider 
the CMB power spectrum, needed to specify the likelihood function of the 
extragalactic temperature anisotropics, as known. 

If a quadrupole spectral distortion in the CMB is definitively measured then 
this formalism can be exploited to constrain the closest stellar encounter. 
There are similiar hopes to constrain such a stellar perturbation by using the 



^■^ Gravitational focusing will induce an additional correlation between the distribu- 
tion of b and Vrei, but we ignore it here since gravitational focusing is not consistent 
with our use of the impulse approximation. 
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extreme orbits of Sedna and analogs that should be discovered in the future. It 
remains to be seen if this goal can actually be accomplished as the associated 
CMB spectral distortion moments suffer from large statistical fluctuations and 
the initial conditions of the Sedna-like objects are unknown. 



6 Conclusions 



We analyzed the effect of a close stellar encounter on the shape of the Oort 
Cloud and the subsequent anisotropic CMB spectral distortions that would 
result. Much about the Oort Cloud is unknown and our calculations in this 
paper are subject to these uncertainities. We have strived to adopt con- 
servative estimates for our mass model of the Oort Cloud and have high- 
lighted the resulting uncertainities. The most likely minimum impact pa- 
rameter over the past 1 Gyr for the present-day local Solar neighborhood is 
b = 2100 AU. Assuming the perturber is a Solar mass star traveling with 
relative velocity 30 km s~^ with respect to the Sun, the perturbation will 
produce frequency-dependent quadrupole moments in the observed CMB of 
C2 = 4.5(3.5) X 10-l^ 6.7(1.1) x 10^^^ 1.1(0.11) x 10"^ at z/ = 30, 353, 545 GHz, 
respectively in the intensity and (temperature) fluctuations. The asphericity 
in the Oort Cloud that leads to our signal will produce spectral distortions 
in the CMB on all angular scales. Our calculation can be extended to the 
octopole {£ = 3), or any higher moment, without much difficulty. Moreover 
the signal will be highly non-Gaussian and will induce correlations between 
the power spectra on different scales. These features may help in the detection 
of the signal. Using data from scales smaller than quadrupole will also help in 
disentangling the signal from the Galactic foreground emission. 

The Planck Satellitt^ will have nine bands of frequency coverage between 
= 30 GHz and 857 GHz, with the u = 353 GHz band the highest with 
reasonable estimated detector noise. Additionally, the recently launched far- 
infrare d satellite Aka. r i, orig inally Astro-F, has frequency coverage down to 



2 THz (jPearson et al.l . I2OOII ) and might be able to search for the signal from 



the closest and hottest Oort Cloud objects. 

Finally, we described a Bayesian method to use observed anisotropic spectral 
distortions to constrain the trajectory of the closest stellar perturber. The 
Bayesian method requires that the direction of the quadrupole, or equiva- 
lently the phase information in the quadrupole moments, does not drastically 
change over timescales comparable to the period of the Oort Cloud objects. 
The associated fluctuations are due to Poisson fluctuations in the number of 
Oort Cloud objects at any position in the outer Solar System. We quantifled 



13 



www.rssd.esa.int/Planck/ 
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these fluctuations and determined that they are important when interpret- 
ing our calculations due to the relatively small number of test particles (10^) 
we used. However, we estimated that the real Oort Cloud has approximately 
10^^ objects and the interpretation of the real quadrupole moments would 
not be affected by these statistical fluctuations. Therefore the direction of the 
quadrupole of CMB spectral distortions, if ever observed, could be used to 
constrain the properties of the closest stellar encounter with the Solar Sys- 
tem. 
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A Angular Moment Definitions 



In this appendix we review the relationships between the various descrip- 
tions of the angular moments of a distribution on a sphere, in particular the 
quadrupole moments which are discussed in this paper. It is advantageous to 
describe the CMB temperature anisotropies in a basis of spherical harmonics 
because of their orthonormality properties on the sphere. A function defined 
over the sphere, A{n), can be expanded in a spherical harmonic basis as 

aem = J £hY;^{h)A{h). (A.l) 

Since the CMB temperature anisotropies are assumed to be described by a 
Gaussian random field the particular values of the aim are not meaningful, only 
the statistical properties of the expansion coefficients carry any information 
about the underlying cosmology. Therefore, it is customary to calculate the 
ensemble average of the a^m covariance matrix 

(fl£imifl£2»Tt2) = ^^i/2'^mi,-m2(~l)'"^C'^i5 (^-2) 

which can be related to an isotropicP^ quantity, the power spectrum, Ci. This 
identity can be used to motivate an observational definition of Ci 



1 



21 + 



_ ^ (-l)™a^„a^_^. (A.3) 



Invariant under rotations of our coordinate system. For example, it is invariant 
with respect to different choices of our reference direction. 
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We will now relate these quantities to the Cartesian multipole moment expan- 
sions utilized in our paper and defined in Eq. (fT^ . Expressing the Cartesian 
directional vectors in spherical coordinates, 



J4> 



fix + ifiy = sin I 
n-. = cos 6. 



Then the five i = 2 spherical harmonics can be expressed as 



(A.4) 
(A.5) 



Y2i{h) 

Y2^iiri) 
Y2-2{n) 



/3 [T^ 



/3 5 



(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A.IO) 



/ 8 V 47r 

From these definitions relationships between the a2m and Qij can be derived. 



/3 5 



^22 



0-21 



020 



02-1 



02-2 




8 V 47r 



i^Qxx Qyy '^^Qxy)i 



2 V 47r ^^^^ ^Qyz); 



3 h_ 
2V 4^ 

V 2 V 47r 



Q ZZ ) 

i^Qxz ~l~ ^Qyz^i 



/3 / 5 

8 V 47r ^^^^ Qyy ~^ '^^Qxy^- 



(A.ll) 
(A.12) 
(A.13) 
(A. 14) 
(A.15) 



The power spectrum can be related to the Cartesian quadrupole moments 
through Eq. ([Qj) 



C2 = -^{Qlx + Qlz + Qly + Qlz + Qlz + QxxQzz)- (A.16) 

The expectation value of this expression should contain both the mean and 
variance of the Qij as terms. However, consistent with our calculation that 
Poisson fiuctuations are not important we will evaluate Eq. (1A.16P with only 
the means of the Qij ■ 
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y ct/Tcmb X 10-*^ FWHM (arcmin) 

30 GHz 2.0 33 

100 GHz 2.5 9.5 

217 GHz 4.8 5.0 

353 GHz 14.7 5.0 

545 GHz 147 5.0 

Table B . 1 " 

Instrument noise parameters adopted for this paper. These values were taken from 

the Planck Science webpage (www.rssd.esa.int/Planck/). 

B Instrument Noise 



Finally, we describe the instrument noise values that were adopted in this 
paper. In tem perature units the instrument noise for the observation of Ci is 
f lKnoxl . Il99,^ 



47rcr2 



where the beam width is related to the FWHM as V8 ln2aB =FWHM, Npi^ 
is the number of observed pixels and a is the noise per pixel in temperature 
units. The temperature fluctuation can be related to an intensity fluctuations 
as follows, 

A/. = ^^|TeMBAT.. (B.2) 
Table [B] lists the noise parameters used in this paper. 
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